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We present a measurement of the partial branching fractions and mass spectra of the exclusive 
radiative penguin processes B — > Kit-kj in the range mx-itn < 1.8 GeV/c 2 . We reconstruct four final 
states: K + ir~ir + ~/, K + 7r _ 7r°7, K^ty~ 7r + 7, and K 3 \w + 'K y, where K% — > 7r + 7r~. Using 232 million 
e + e~ — > BB events recorded by the BABAR experiment at the PEP-II asymmetric-energy storage 
ring, we measure the branching fractions B(B + — ► K + ir~ ir +r y) = (2.95±0.13 (stat.)±0.20 (syst.)) x 
10" 5 , B(B° -> K+ir-Tr -/) = (4.07±0.22 (stat.)±0.31 (syst.)) x 10~ 5 , B(B° -> K°-K + ir~~/) = (1.85± 
0.21 (stat.)±0.12 (syst.)) x 10~ 5 , and B(B+ -f K"-k + -k ^) = (4.56±0.42 (stat.)±0.31 (syst.)) x 10~ 5 . 



PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 

In the standard model (SM) the radiative penguin de- 
cay B — > X s 7, where X s is a hadronic system with 
unit strangeness, proceeds via weak-interaction loop dia- 
grams. New physics, beyond the SM, may also contribute 
to the loop amplitude, and lead to differences from the 
SM. This possibility has been pursued in inclusive mea- 
surements, which are theoretically clean but experimen- 
tally challenging, and in exclusive measurements, such 
as B — > Kix^j. We report measurements of the branching 
fractions and mass spectra for the decays B — > Kmrj 
in four channels. SM predictions of the rates and res- 
onance structure of these decays have large uncertain- 
ties The K + ir + ir~"f and A'°7r + 7r~7 decay channels 
have previously been observed ■ Throughout this Let- 
ter, stated decays include charge conjugate modes. 

The decays B — > Kit + ir°"j, which have not previ- 
ously been observed, are of particular interest because 
these three-body hadronic states permit the measure- 
ment, given sufficient statistics, of the photon polar- 
ization The polarization measurement depends on 
the interference between processes such as (A"7r + )7r°7 
and (_ftT7r°)7r + 7, where () indicates resonant substructure. 
This measurement may be compared with the SM pre- 
diction of nearly complete left-handed polarization. 

We use a sample of (232±1.5) x 10 6 BB pairs in a 210.9 
fb _1 dataset collected at the T(45) resonance with the 
BABAR detector at the PEP-II asymmetric-energy e + e~ 
collider. For background studies, we also use a 21.7 fb _1 
sample collected below the BB threshold. The measure- 
ment procedure was designed using simulated signal and 
background events, data in sideband kinematic regions, 
and reconstructed B — » Dtt + , D — > Ktttt decays. Only 
after we established the selection and fit procedures did 
we examine signal candidates in the data sample. 

A description of the detector exists elsewhere For 
this measurement, the most important detector elements 
are the five-layer silicon microstrip tracking detector 
(SVT) and the forty-layer drift chamber (DCH), situ- 
ated in a 1.5 T solenoidal magnetic field, which measure 
charged particle momenta; the CsI(Tl) electromagnetic 
calorimeter (EMC), which measures the energies and di- 
rections of the photons; and the detector of internally 
reflected Cherenkov light (DIRC). The DIRC response 



and energy loss (dE/dx) measured in the SVT and DCH 
are used to identify charged kaons and pions. 

We reconstruct the photon candidate in the Knit"/ de- 
cay from an EMC shower not associated with a charged 
track. The photon must be in the fiducial region of the 
EMC, have a shower- profile consistent with a single pho- 
ton, and be well-separated from other showers. To re- 
move photons from 7r° (rf) decays, we combine the can- 
didate with other photons having energies of at least 50 
(250)MeV/c 2 , and reject it if the invariant mass of any 
combination is within 25 (40) MeV/c 2 of the ir° (77) mass. 

We select if ± and tt^ candidates from charged tracks 
consistent with a kaon or pion mass hypothesis in the 
DIRC and in the dE/dx in the SVT and DCH. We re- 
construct K® candidates from pairs of oppositely-charged 
tracks, and determine the decay vertex with a fit. We re- 
quire that the invariant mass falls within 1 1 MeV/c 2 of 
the K® mass; that the distance between the B decay ver- 
tex and the K ° vertex exceeds 5 times the uncertainty 
on the distance; and that the angle between the K® tra- 
jectory and its momentum is less than 100 mrad. We re- 
construct 7T° candidates from pairs of EMC showers each 
with energy >50MeV. We require the invariant mass to 
be within 16 MeV/c 2 of the ir° mass, and that the energy 
of each pair in the T(4S) center of mass (CM) frame ex- 
ceeds 450 MeV; this last selection is about 83% efficient. 

The dominant source of background is continuum pro- 
duction of light quark-antiquark pairs, in which a high- 
energy photon typically is produced either by initial state 
radiation, or from the decay of a ir° or 77 in which one 
photon is not detected. To reject these backgrounds, we 
construct a Fisher discriminant [f| from the polar angle of 
the B candidate in the CM frame, the angle between the 
thrust axis of the B and the thrust axis of the remaining 
charged and neutral particles, and the ratio of the second 
to zeroth angular moments of the remaining charged and 
neutral particles around the thrust axis of the B. We op- 
timize the coefficients independently in each channel to 
discriminate between simulated signal and continuum. 

We perform a geometric fit to the reconstructed B can- 
didate, with production vertex constrained to the nomi- 
nal beam spot, rejecting the candidate if the final state is 
inconsistent with decay from a single vertex. We define 
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AA* = E* B - A* cam and m ES = V^bLn - P$, where 
E B and p*j are the CM energy and momentum of the 
B candidate, and AjJ eam is the CM energy of each beam. 
We require m ES > 5.2GeV/c 2 and |AA*| < 0.15 GeV. 
We also require that the invariant mass of the Kirn sys- 
tem, rnKn-K, fall below 1.8 GeV/c 2 ; this eliminates much of 
a rising continuum background with very little expected 
signal loss. It also removes Ktttt combinations from D 
decays, in B — > Dir° and B Dr\ where the ir° or r\ are 
mis-reconstructed as a photon. In an event in which we 
reconstruct multiple candidates in one channel that pass 
the selection requirements (occurring in 11-27% of se- 
lected signal events, depending on the channel), we keep 
the candidate with the largest vertex probability (with 
the best 7T° mass in the A°7r + 7r°7 channel; with the ir° 
mass as a tie breaker in the K + ir~ir°~f channel) and reject 
the others. Candidates reconstructed in different chan- 
nels are allowed in the same event. The dependence of 
the efficiency of our selection requirements on intermedi- 
ate resonance and on ms;„ has been checked and found 
to be small; systematic uncertainties are discussed later. 

The dominant backgrounds from BB events after the 
selection criteria have been applied are b —* 57 processes. 
We categorize these backgrounds: (i) "crossfeed" from 
mis-reconstructed Kttttj decays, such as by choosing in- 
correctly a particle from the other B; (ii) B — > Kttj 
decays that combine with a track from the other B to 
form a Kimj candidate; and (Hi) backgrounds from all 
other b — ► sj decays. A crossfeed candidate may be re- 
constructed in the same decay channel in which it is pro- 
duced, or in a different channel, and can also be produced 
in a B — > Kmr^/ decay that is not used in this analysis 
(such as B + — > A' + 7r°7r°7). We model our signal as well 
as crossfeed backgrounds with simulated B — > A'^7 de- 
cays, where Kx is any of the five lowest-lying J > 
kaon resonances above the A* (892). We study back- 
grounds from A7T7 using simulated B — > A* (892)7 and 
B — * A| (1430)7 decays. We study backgrounds from 
other b — > 57 decays using an inclusive simulation ac- 
cording to the model of Kagan and Neubert @ with 
mj) = 4.8 GeV/ c 2 , tuned to match multiplicity distribu- 
tions measured in inclusive b — > sj decays Q ■ The largest 
final background contributions from b — > sj processes are 
crossfeed backgrounds, for which we obtain yields rang- 
ing from 55% to 95% of the signal yields. 

We estimate other sources of background candidates 
from B decays other than b — * S7 processes by simu- 
lating generic B decays. We pay special attention to B 
decays with Kitttit and Kirirq final states; if the 7r° or 
r\ decays asymmetrically and we don't detect the lower- 
energy photon, the kinematic properties of the resulting 
B candidate may resemble a signal candidate. We study 
these decays using high-statistics simulated samples, and 
look for signal candidates that are reconstructed from a 
single B — > Kit-kit or B — > Kirirr] decay. We expect to 
reconstruct fewer than two such candidates per channel. 



We perform a maximum likelihood fit to the joint 
toes-AA* distribution of our selected candidates. We fit 
all four channels simultaneously to account for crossfeed 
backgrounds between channels. The likelihood function 
contains terms for correctly reconstructed signal candi- 
dates, crossfeed background candidates between all 16 
combinations of the production and reconstruction chan- 
nels, backgrounds from B — ► A-7T7 and from other b — > sj 
decays, and backgrounds from continuum events. We 
have determined from simulations that the dominant con- 
tinuum background component adequately accounts for 
combinatoric backgrounds from other BB decays, which 
do not show strong peaks in toes and A A* . 

The likelihood function for a candidate reconstructed 
in decay channel i with kinematic variables y = 
(m E s, AA*) is given by, 

fun = n bb U^jUy) + &<£f?(v)j 

+ <fi(y)+nifi(y) , 

where N BB is the number of BB pairs in our dataset; B % 
is the branching fraction for decay channel i; t\ and j\ are 
the efficiency and probability density function (PDF) for 
correctly reconstructed signal candidates in decay chan- 
nel i; and fl l are the efficiency and PDF for crossfeed 
background candidates produced in channel j and recon- 
structed in channel i; n l c and /* are the yield and PDF 
for backgrounds from continuum and generic BB decay 
events in channel i; and n\ and are the yield and PDF 
for backgrounds from other b — > 37 processes in channel 
i. We further parameterize the likelihood function by the 
four data-taking runs during which data were collected, 
accounting for slight changes in experimental conditions. 

The branching fractions B l , yields n*, and shape pa- 
rameters of fl are varied in the fit; other efficiencies, 
yields, and PDF shapes are fixed from simulation stud- 
ies. We parameterize /] as the product of Crystal Ball 
functions § of toes an d of AA*, /* as the product of 
a Crystal Ball function of toes and a linear function of 
AA*, and /* as the product of an Argus function [§] 
of Toes and an exponential function of AA*. We use a 
binned parameterization for /£. As the signal and cross- 
feed terms are both scaled by the parameters B l , the 
crossfeed background yields vary with the signal branch- 
ing fractions, and we measure the branching fractions 
from yields of both signal and crossfeed candidates. 

Tablc[T]shows the fit results. Projections in toes, along 
with the fit results, are displayed in Fig. [1] The fit proba- 
bility (P- value) is evaluated with a likelihood ratio statis- 
tic [10], assuming Poisson-distributed bin contents, to be 
10%. The distribution of the test statistic under the null 
hypothesis is evaluated by simulation. 

Figure [5] shows background-subtracted ttik-kt: mass 
spectra. Background subtraction is achieved using the 
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FIG. 1: Distributions of tties (points). Projected tties distri- 
butions from the fit are shown as cumulative curves: contin- 
uum and generic BB component (dashed), b —* S7 component 
(dotted, includes cross- feed), and signal (solid). The small os- 
cillation in the dotted and solid curves is due to the use of 
binned distributions to model the b —* S7 component. 



results of the fits to calculate event-by-event weights to 
extract the signal component We present branch- 

ing fractions in bins of rriKTrn, which are largely model- 
independent, instead of extracting B — > Kx^f branching 
fractions for specific Kx resonances. Disentangling the 
resonance structure requires careful modeling of ampli- 
tudes and relative phases of interfering processes, includ- 
ing in the decays of the Kx resonances, not all of which 
are well measured. A partial wave analysis to extract the 
resonance structure and measure the photon polarization 
should be possible with future datasets. 



TABLE I: Results of the fit for B — > Knnj, for m K ww < 
1.8 GeV/c 2 . The first error is statistical, the second system- 
atic. The yields do not include the channel crossfeeds, which 
are included in the fit to obtain the branching fractions. 

Channel Yield Branching Fraction (10 ) 



K + w~ir + -( 899 ±38 



K+T 
K°7V + 



7T°7 



572 ± 31 



7 176 ±20 
7T°7 164 ±15 



2.95 ± 0.13 ± 0.20 
4.07 ± 0.22 ± 0.31 
1.85 ± 0.21 ± 0.12 
4.56 ± 0.42 ± 0.31 



FIG. 2: Background-subtracted vtik^t: spectra. The branch- 
ing fraction in each bin is computed from the weighted event 
yield. Error bars show statistical uncertainties; the systematic 
uncertainties due to b — -> S7 model assumptions are small. 



We validate the procedure for extracting branching 
fractions and mx-mr distributions using fits to simulated 
samples. We verify that the branching fractions and mass 
spectra obtained from these toy fits reproduce on average 
the simulation inputs. We use the same procedure to ex- 
tract the rnKir-n distributions for continuum and generic 
BB backgrounds and for backgrounds from b — > sj de- 
cays; these are consistent with the expected distributions. 

Systematic uncertainties arise from various sources, 
shown in Table [Til The largest sources are: (i) The 
T(4S) branching fractions to B + B~ and B°B° are each 
assumed to be 0.5. We assign a 2.6% systematic uncer- 
tainty to this, based on current information [l2l ]. (ii) The 
uncertainty on the photon selection efficiency determined 
from simulated events is estimated to be 2.7%. (Hi) From 
studies of B — > Dir , D — > Kirir events, we assign an un- 
certainty of 4.2% to the charged kaon identification effi- 
ciency, (iv) The uncertainty of the ir° selection efficiency 
is estimated at 3.0%. (v) There is considerable uncer- 
tainty in the models we use to estimate backgrounds, 
including cross-feed dependence, from b — ► processes. 
We estimate the effect of this uncertainty on both the 
branching fractions and mass spectra by simulating these 
backgrounds with substantially different models. The 
largest effect is in the K®ir~ir + j channel, where the un- 
certainty is 4.0%. (vi) We measure a shift in the beam 
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energy in B — > Dn ± decays, on average 0.6 MeV; we es- 
timate the effect of this on our fits, (vii) We estimate 
bias in the fit due to uncertain parameterization of the 
signal and background PDFs. The largest effect is in the 
if°7r + 7r°7 channel, where the uncertainty is 3.5%. 

We have measured branching fractions for B — > Ktttt^ 
in four decay channels for mKm < 1.8GeV/c 2 . The 
Ktt + t:~ channels are consistent with the previous mea- 
surement [2J. We present first observations of decays in 
the _KT7r + 7r°7 channels that are important to measuring 
the photon polarization. The branching fractions are rel- 
atively large in the context of B — > X s j decays, provid- 
ing encouragement that a polarization measurement may 
be possible with future datasets. Mass spectra for the 
Kinr system are also presented. We observe an enhance- 
ment near 1.3 GeV/c 2 and substantial branching fractions 
at higher masses. Untangling the resonant contributions 
presents a challenge for the polarization measurement. 
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TABLE II: Estimated systematic uncertainties in the branch- 
ing fractions, in percent, by source and decay channel. 



Source 








T AI + (1 

K gTT 7T 


BB count 


1.1 


1.1 


1.1 


1.1 


T(45) BF 


2.6 


2.6 


2.6 


2.6 


Efficiencies: 










Photon selection 


2.7 


2.7 


2.7 


2.7 


7T° and 77 veto 


1.0 


1.0 


1.0 


1.0 


Tracking 


1.4 


1.1 


1.1 


0.8 


Kaon selection 


4.2 


4.2 


1.6 


1.6 


7v^~ selection 


1.4 


1.0 


1.4 


1.0 


7T° selection 




3.0 




3.0 


Fisher cut 


1.0 


1.0 


1.0 


1.0 


Vertex probability 


0.7 


0.7 


0.7 


0.7 


MC statistics 


0.5 


0.7 


0.8 


1.1 


Backgrounds: 










b — > S7 model 


1.4 


1.0 


4.0 


1.3 


B -> KmvK°/rj 


0.2 


0.1 


0.0 


0.6 


Beam energy shift 


1.0 


0.5 


1.6 


0.6 


PDF shape 


0.1 


2.9 


0.9 


0.2 


Fit bias 


1.6 


1.3 


1.3 


3.5 


Total 


6.7 


7.6 


6.7 


6.8 



We present a measurement of the partial branching fractions and mass 
spectra of the exclusive radiative penguin processes B — > Kmry in the range 
rriKixTx < 1.8GeV/c 2 . We reconstruct four final states: K + tt~it + j, K + Tr~n a j, 
K®ir~Tr + j, and K®ir + ir~j, where K® — > n + Tr~ . Using 232 million e + e~ — > BB 
events recorded by the BaBar experiment at the PEP-II asymmetric-energy stor- 
age ring, we measure the branching fractions B(B + — > K + n~ tt + 7) = (2.95 ± 
0.13 (stat.) ± 0.20 (syst.)) x HT 5 , B{B° -> if+Tr-Tr^) = (4.07 ± 0.22 (stat.) ± 
0.31 (syst.)) x 10~ 5 , B(B° -> K°ir + %-'y) = (1.85 ± 0.21 (stat.) ± 0.12 (syst.)) x 
10- 5 , and B{B+ -> if°7r+7r°7) = (4.56 ± 0.42 (stat.) ± 0.31 (syst.)) x 10~ 5 . 
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